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Gaussian bipartite states are basic tools for the realization of quantum infor- 
mation protocols with continuous variables. Their complete characterization 
is obtained by the reconstruction of the corresponding covariance matrix. Here 
we describe in details and experimentally demonstrate a robust and reliable 
method to fully characterize bipartite optical Gaussian states by means of a 
single homodyne detector. We have successfully applied our method to the 
bipartite states generated by a sub-threshold type-II optical parametric os- 
cillator which produces a pair of thermal cross-polarized entangled CW fre- 
quency degenerate beams. The method provide a reliable reconstruction of the 
covariance matrix and allows to retrieve all the physical information about the 
state under investigation. These includes observable quantities, as energy and 
squeezing, as well as non observable ones as purity, entropy and entanglement. 
Our procedure also includes advanced tests for Gaussianity of the state and, 
overall, represents a powerful tool to study bipartite Gaussian state from the 
generation stage to the detection one. 
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1. Introduction 

The quantum characterization of physical systems has a fundamental interest in its own 
and represents a basic tool for the design of quantum protocols for information processing in 
realistic conditions. In particular, the full experimental reconstruction, at the quantum level, 
of optical systems opens the way not only to high fidelity encoding/transmission/decoding 
of information, but also to the faithful description of real communication channels and to 
precise tests of the foundations of quantum mechanics [1-4]. 

Among the systems of interest for quantum information processing we focus on the class 
of bipartite optical states generated by parametric processes in nonlinear crystals. These are 
Gaussian states and play a crucial role in quantum information processing with continuous 
variables [5-8]. Indeed, using single- and two- mode Gaussian states, linear optical circuits 
and Gaussian operations, like homodyne detection, several quantum information protocols 
have been implemented, including teleportation, dense coding and quantum cloning [9]. In 
particular, Gaussian entangled states have been successfully generated in the laboratories by 
type-II optical parametric oscillators (OPO) below threshold [10-14]. In these OPO systems 
the parametric process underlying the dynamics is well described, at least not too close to 
the threshold, by bilinear Hamiltonian, thus the output states are Gaussian and they are 
completely characterized by the first and second moments of their quadratures, i.e. the 
covariance matrix. 

In this paper we address characterization of bipartite Gaussian states and review in 
details a scheme to fully reconstruct the Gaussian output from an OPO below threshold, 
which has been proposed in the recent years [15, 16] successfully implemented experimentally 
[17]. In the present contibution we give a more accurate description of the experiment and 
data analisys and, in particular, we pay attention to advanced Gaussianity test beyond the 
simple check of Kurtosis. Our method relies on a single homodyne detector: it provides the 
full reconstruction of the covariance matrix (CM) by exploiting the possibility of optically 
combining the two frequency degenerate OPO signal and idler beams and then measuring 
suitable quadratures on the obtained auxiliary modes. Once the CM is obtained one may 
retrieve all the quantities of interest on the state under investigation, e.g. energy and 
squeezing, including those not corresponding to any observable quantity like purity, entropy, 
entanglement, and mutual information. Quantum properties are discussed in view of the 
possible use of these states in quantum communication protocols. In particular, we address 
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the dependence of mutual information as a function of the bi-partite system total energy. 
Of course, a bipartite state is fully characterized by its covariance matrix if and only if it is 
a truly Gaussian one. Usually one assumes that the state to be processed has a Gaussian 
character because the interaction Hamiltonians are approximated by bilinear ones and this 
is often an excellent approximation [18]. In turn, the resulting evolution corresponds a 
Gaussian operations. On the other hand, it is known that nonGaussian dynamics may 
occur when the OPO approaches the threshold [19, 20] and when phase diffusion [21, 22] 
is present during the propagation and/or the detection stages. Therefore, in order to avoid 
any possible experimental issue [23-26], a preliminary check on the Gaussian character of 
the signal is crucial to ensure that the actual measured CM fully characterizes the quantum 
state. For the first time, in this paper, CM data analisys includes advanced statistical tests 
to assess Gaussianity [19, 27] of the state. The complete characterization strategy represents 
a powerful tool to study bipartite Gaussians state from the generation stage to the detection 
one. 

The paper is structured as follows. In section 2 we introduce the formalism used through- 
out the paper and, in particular, we review two-mode Gaussian states and their covariance 
matrix as well as the relations among the CM elements and some physical quantities of 
interest, such as the purity, the entropy and the entanglement. The method to reconstruct 
the CM is described in section 3, while section 4 is devoted to the details of our experimental 
implementation. The analysis of the data and the results are discussed in details in Sections 
5, 6 and 7. In particular, test of Gaussianity are illustrated in Section 5 and results from 
full quantum tomography in Section 6. Section 8 closes the paper with some concluding 
remarks. 

2. Two-mode Gaussian states 

A n-mode state g of a bosonic system is Gaussian if its characteristic function = 
Tr[^D(A)] has a Gaussian form, D(X) = ®J! =1 Dfe(Afe) being the n-mode displacement op- 
erator with A = (Ai, . . . , A n ), Afc e C, and D k (X k ) = exp{X k a\. — A^a^} denoting single-mode 
displacement operators [7]. Gaussian states are completely characterized by the first and 
second statistical moments of the quadrature field operators, i.e. by the vector of mean val- 
ues and by the covariance matrix (CM). Since in this paper we focus on two- mode Gaussian 
states of the radiation field, in this section we review the suitable formalism to describe the 
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system. We also assume, since this is the case in our experimental implementation, that 
the mean values of quadratures are zero. Upon introducing the vector of canonical oper- 
ators R = {x a ,y a ,Xb,yb), in terms of the mode operators a k , k = a,b, x k = + a k ), 
jjk = — Ofc) the CM cr of a bipartite state is the real symmetric definite positive block 
matrix: 



A 



C 



, . (!) 

C 

with /Jhk = U{R k ,R h }) - (R k ){R h ) being {f,g} = fg + gf. Matrices A, B and C are 2 x 2 
real matrices, representing respectively the autocorrelation matrices of mode a and b and 
their mutual correlation matrix. It can be observed that each block A, B and C can be 
written as the sum of two matrices, one containing the product of mean values of (R k ) and 
the other contain the mean value of products of operators {{Rk, Rh}) ■ 

Once the CM is known, all the properties of g may be described and retrieved. As for 
example, the positivity of g, besides positivity of the CM itself, impose the constraint 

<r + l -n > o, (2) 

where f2 = u> © u> is the two- mode symplectic matrix, given in terms of u> = adiag[l, —1]. 
Inequality (2) is equivalent to the Heisenberg uncertainty principle and to positivity, and 
ensures that cr is a bona fide CM. 

A relevant result concerning the actual expression of a CM is that for any two-mode CM 
cr, there exists a (Gaussian) local symplectic operation S — Si © S 2 that brings cr in its 
standard form, namely [28, 29]: 



S t ctS 



A 


C 


C T 


B 



(3) 



where A = diag[n, n], B = diag[m, m], C = diag[ci, C2], with n, m, c\ and C2 determined by 
the four local symplectic invariants l\ = det(A) = n 2 , I2 = det(B) = m 2 , I3 = det(C) = 
C1C2, h = det(cr) = {nm — cf) {nm — c|). If n — m, the matrix is called symmetric and 
represents a symmetric bipartite state where the energy is equally distributed between the 
two modes. 

By using the symplectic invariants the uncertainty relation (2) can be expressed as: 

J 1 + / 2 + 2/ 3 <4/ 4 + i (4) 



It is useful to introduce the symplectic eigenvalues, denoted by d± with cL < d+, which in 
terms of symplectic invariants read as follows [30] 



d± = ^ AWiffi , (5) 
where A(<r) = Ji + J 2 + 2/ 3 . In this way, the inequality (2) re-writes as: 

d- > 1/2. (6) 

A real symmetric definite positive matrix satisfying gL > 1/2 corresponds to a proper CM 
i.e. describes a physical state. 
A. Purity and entropies 

The purity of the two- mode Gaussian state g, may be expressed as a function of the CM (1) 
as follows [31]: 

/z^ = Tr[£ 2 ] = (16/ 4 n. (7) 

Another quantity, characterizing the degree of mixedness of g, is the von Neumann entropy 
S (g) = — Tr (g log g). If the state is pure the entropy is zero (S — 0), otherwise it is positive 
(S > 0) and for two-mode Gaussian states it may be written as [7, 30]: S(g) = S(cr) = 
f{d + ) + f(d-) where the symplectic eigenvalues d± are given in (5) the function f(x) = 
(x + 1/2) \og(x + 1/2) — (x — 1/2) \og(x — 1/2). It is useful to recall that for a single mode 
Gaussian state the von Neumann entropy is a function of the purity alone [32]: 

whereas for a two-mode state all the four symplectic invariants are involved. 

For a two-mode state g it is of interest to assess how much information about g one can 
obtain by addressing the single parties. This is of course related to the correlation between 
the two modes and can be quantified by means of the quantum mutual information or the 
conditional entropies [33]. Given a two-mode state g the quantum mutual information 1(g) 
is defined starting from the von Neumann entropies as: 

I(g) = S(g l ) + S(g 2 )-S(g), 

where g^ = Tr h (g), with k,h = 1,2 and h ^ k, are the partial traces, i.e. the density 
matrices of mode k, as obtained tracing over the other mode. / (g) can be easily expressed 
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in terms of the blocks of cr and its symplectic eigenvalues. One has 



1(a) = f (y/h) + f (v 7 ^) - f(d+) - f(d-), 



(9) 



where f(x) is reported above. The conditional entropies are defined accordingly as [33]: 



If 5'(1|2) > (or 5(211) > 0), the conditional entropy gives the amount of information 
that the party 1 (2) should send to the party 2 (1) in order to allow for the full knowledge 
of the overall state g. If ^(l^) < (or S(2\l) < 0), the party 1 (2) does not need to 
send any information to the other and, in addition, they gain — S'(l|2) or — 5'(2|1) bits of 
entanglement, respectively. This has been proved for the case of discrete variable quantum 
systems [34] and conjectured [35] for infinite dimensional ones. 
B. Entanglement 

A two-mode quantum state p is separable if and only if it can be expressed in the following 

form: p = J2kPk {pt^ ® Pk^i wrt h Pk > 0, J2kPk = 1 an d pi^ an d pf 1 are single-mode 
density matrices of the two modes a and b, respectively. Viceversa if the state is not sep- 
arable, it is entangled. A general solution to the problem of separability for mixed state 
has not been found yet. For two- mode Gaussian states there exist necessary and sufficient 
conditions to assess whether a given state is entangled or not. In particular, there are two 
equivalent criteria, usually referred to as Duan criterion and Peres- Horodecki- Simon crite- 
rion, which found an explicit form in term of the CM elements. The criteria provide a test 
for entanglement, whereas to assess quantitatively the entanglement content of a state one 
may use the logarithmic negativity or the negativity of the conditional entropies, as we see 
below. 

1. Duan criterion 

This criterion [29] is based on the evaluation of the sum of the variances associated to 
a pair of EPR-like operators defined on the two different subsystems. For any separable 
continuous variable state, the total variance is bounded by twice the uncertainty product. 
For entangled states this bound can be exceeded and the violation provides a necessary 
and sufficient condition for entanglement. The criterion leads to an inequality that can be 



S(2\1) = S{q)-S{q 1 ). 



S(l\2) = S(q) - Sin*) 



(10) 
(11) 
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expressed in terms of standard form CM elements: 



Tft 1 

/3 D = na 2 + — - |ci| - \c 2 \ < a 2 + — , (12) 
a 2 a 2 



with a 2 = y ^zj- A separable state will not satisfy the above inequality. The criterion 
raises from the fact that for an entangled state it is possible to gain information on one of 
the subsystems suitably measuring the other one. 

2. Peres-Horodecki-Simon criterion (PHS) 

Also PHS criterion establishes a necessary and sufficient condition for separability of bipar- 
tite Gaussian states [28]. Given the CM cr, the corresponding two- mode Gaussian state is 
not separable iff: 

a- + ^n<o, (13) 

where A = diag[l, 1, 1, —1] and cr = AcrA is the CM associated with the partially trans- 
posed density matrix. Thanks to the symplectic invariants {ii, I 2 , I3, I4} the inequality 
(13) can be written in a form that resembles the uncertainty relation (4) [7]: 

I 1 + I 2 + 2\I 3 \>4h + ^, (14) 

or, in terms of standard form CM elements, as: 

1 



or simply as: 



where : 



n 2 + m 2 + 2 \cic 2 \ — 4 {nm — cfj (nm — c^) < -, (15) 



J_ < 1/2, (16) 



d ± = 



A(tr) ± \ A(<r) 2 - 4J 4 
\ 2 » ( 17 ) 

are the symplectic eigenvalues of & and A(c) = l\ + I 2 — 2/3. Therefore, iff oL < 1/2 the 
Gaussian state under investigation is entangled. 

For an entangled state a quantitative measure of entanglement can be given on the 
observation that the larger is the violation gL < 1/2 the stronger the entanglement, or 
more properly, the stronger the resilience of entanglement to noise [36-39] . The logarithmic 
negativity for a two-mode Gaussian state, is given by [40] (remind that gL > 0) 

E{&) = max j(D,-log2cLj , (18) 



and it is a simple increasing monotone function of the minimum symplectic eigenvalue cL (for 
gL < 1/2): it thus represents a good candidate for evaluating entanglement in a quantitative 
way. It is worth to note that also the negativity of the conditional entropies (10) and (11) 
is a sufficient condition for entanglement [41]. 

C. EPR correlations 

This way of assessing quantum correlations between two modes is named after the analogy 
with the EPR correlation defined for a system undergone to a quantum nondemolition 
measurement (QND)[13]. Let us consider two subsystems a and b, QND establishes, in 
principle, that measurement performed on subsystem b, does not affect system a. This 
criterion is equivalent to state that the conditional variance V a \ b of a quadrature of beam 
a, knowing beam b, takes a value smaller than the variance a would have on its own. The 
conditional variance can be expressed in terms of the unconditional variance V a of subsystem 
a (i.e. the variance that the same quantity has in the subsystem a space) and normalized 
correlation C a b between the two [14, 42]: V a \ b = V a (l — C^ b ) , an analogous relation holds 
for Vb\ a . A bipartite state is said to be EPR correlated if it verifies the following inequality: 

V a \ b V b \ a < 1/4, (19) 

that can be rewritten in terms of standard form CM elements as follows: (5e = 
nm ( 1 — ^- ) ( 1 — ^- ) < 1/4. If the inequality is satisfied in the system described by 

y nm I y nm I I ~l j j j 

the CM the information on a (b) extracted from a measurement on b (a) is sufficient for 
knowing its state with a precision better than the limit given by the variance for a coherent 
state. In turn, EPR correlations are stronger than entanglement [10, 24, 43], i.e. all EPR 
states are entangled whereas the converse is not true and there are entangled states violating 
Ineq. (19). 

3. Covariance matrix reconstruction 

In this section we describe in some detail the method we have implemented to experi- 
mentally reconstruct the CM given in Eq. (1) and, thus, to fully characterize a bipartite 
Gaussian state. As expected, each autocorrelation block, A or B, is retrieved by mea- 
suring only the single-mode quadratures of the concerned mode a or b. Diagonal terms 
of A correspond to the variances of x a and y a and are directly available at the output of 
the homodyne detection. Off-diagonal terms are instead obtained by measuring the two 
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additional quadratures z a = ^ (x a + y a ), t a = (x a — y a ), and exploiting the relation 
o"i2 = 021 = \{( z 1) ~ (O) ~~ ( x a)(y a ) [15]. The block B is reconstructed in the same way 
from the quadratures of b. Elements of block C involve the products of quadrature of modes 
a and b and cannot be obtained by measuring individually the two modes. Instead they are 
obtained by homodyning the auxiliary modes c = -^(a + b), d = ^(a — 5), e = -^{ia + b), 
and / = -j=(ia — b) and by making use of the following relations: 

*13 = \{{ x l) - { X D) - (X a ){x b ), 
<?U = \{{VI) ~ (Vj)) - (Xa){Vb), 
^23 = \(( X )) - ( X D) - (ya)( x b), 

°2A = ^{(y 2 c ) - (y 2 d )) - (y a )(y b ) 

It is worth to note that since (x 2 f ) = (xf) + (y a ) 2 — (x 2 e ) and (yj) = (x 2 a ) + (y 2 ) — (y 2 ), the 
measurement of the /-quadratures is not mandatory. 

As we will see in the following, our experimental setup allows one to mix the modes a 
and b, say the signal and idler, thank to the polarization systems at the OPO output. At 
the same time, the quadratures x = xq, y = x n /2, z = x^/i and t = x~ n /4 required for the 
entanglement measurement and for the reconstruction of the CM can be easily and reliably 
reconstructed from the pattern function tomography applied to data collected in a 2n scan 
of the homodyne detector. 

4. Experimentals 

The experimental setup is schematically depicted in Fig. 1. It is based on a CW internally 
frequency doubled Nd:YAG laser (Innolight Diabolo) whose outputs @532nm and @l064nm 
are respectively used as the pump for a non degenerate optical parametric oscillator (OPO) 
and the local oscillator (LO) for the homodyne detector. The OPO is set to work below the 
oscillation threshold and it provides at its output two entangled thermal states (the signal, 
a and the idler b): aim of the work is indeed to measure the covariance matrix of these two 
beams. 

The OPO is based on an a-cut periodically poled KTP non linear crystal (PPKTP, 
Raicol Crystals Ltd. on custom design) which allows for implementing a type II phase 
matching with frequency degenerate and cross polarised signal and idler beams, for a crystal 
temperature of ~ 53°C. The transmittivity of the cavity output mirror, T out , is chosen in 
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order to guarantee, together with crystal losses (k) and other losses mechanisms (T in ), an 
output coupling parameter i] out = T out /{T in + k) @1064 nm of ~ 0.73, corresponding to an 
experimental linewidth of 16 MHz @1064 nm. In order to obtain a low oscillation threshold, 
OPO cavity geometry is set to warrant simultaneous resonance on the pump, the signal 
and the idler: pump resonance is guaranteed by servo-assisting the OPO cavity with a 
Drever Pound Hall system [44] , while the resonance of other beams is induced by exploiting 
the natural birefringence of the KTP to tune the optical path of each beam inside the 
cavity, through a fine control of the crystal temperature and tilt [45]. The OPO is equipped 
with an handmade control system able to stabilise the non-linear crystal temperature up 
to O.lmK. Measured oscillation threshold is 50 mW; during the acquisition the system has 
been operated below threshold at ^60% of the threshold power. 

The signal and idler modes are then sent to the covariance matrix measurement set- 
up: this consists in a preliminary polarisation system, that allows choosing the beam to 
be detected and a standard homodyne detector. The polarisation system is made of an 
half- wave plate (A/2) followed by a polarising beam splitter (PBS); the different wave-plate 
orientations allow choosing the beam to be transmitted by the PBS: the signal (a), the idler 
(b) or their combinations c and d. The other auxiliary modes e and / may be obtained by 
inserting before the PBS an additional quarter wave plate (A/4) [15]. Acquisition times are 
considerably short thank to pc-driven mechanical actuators that allow setting the A/2 and 
A/4 positions in a fast and well calibrated manner. 

Once a beam is selected, it goes to a homodyne detector put downstream the PBS. This 
exploits, as local oscillator, the laser output @1064 nm, previously filtered and adjusted 
to match the geometrical properties of the OPO output: a typical interferometer visibility 
is 0.98. The LO oscillator phase 9 is spanned thanks to a piezo-mounted mirror, linearly 
driven by a ramp generator which is, in turn, adjusted to obtain a 2n variation in 200 ms. 
The homodyne photodiodes (PDs, model Epitaxx ETX300 ) have both nominal quantum 
efficiencies of «0.91 and each is matched to a low-noise trans-impedence AC (> few KHz) 
amplifier. The difference photocurrent is eventually further amplified by a low noise high 
gain amplifier (MITEQ AU 1U2). 

In order to avoid low frequency noise, the photocurrent is demodulated with a sinusoidal 
signal of frequency fi=3 MHz and low-pass filtered (5=300 KHz), to be sent to a PCI 
acquisition board (Gage 14100) that samples it with a frequency of 10 6 pts/run, with 14-bit 
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resolution. The total electronic noise power of the acquisition chain is 16 dBm below the 
shot noise level, corresponding to the a SNR ^40. 

5. Gaussianity tests 

Since the covariance matrix contains the full information only for Gaussian states, a prelim- 
inary check on the Gaussian hypothesis is necessary in order to validate the entire approach. 
At first, in order to asses the Gaussianity of our data set we have evaluated the Kurtosis 
excess (or Fisher's index) is calculated. Then, once Gaussianity is proved, a more sophis- 
ticate test is used to check the statistical quality of the collected data. In particular, the 
Shapiro-Wilk [46] test checks whether the collected data come from a truly random nor- 
mal distribution, i.e. whether or not the data ensemble is a faithful replica of a Gaussian 
statistical population. We underline the importance of Gaussianity tests, which is usually 
assumed rather than actually verified experimentally on the basis of analysis of OPO data. 

The Kurtosis is the distribution fourth order moment, and can be seen as a sort of 
"peakedness" measurement of a random probability distribution. Compared to the Gaussian 
value of 3<t 2 (where a is the standard deviation) the Kurtosis-excess 7 is defined as 

where x is the mean of the data and Pi is the probability of the i-th outcome. A 7 = 
distribution is Gaussian. As a matter of fact 7 gives an immediate check on the Gaussian- 
ity of the data ensemble, whereas it cannot say anything about accidental (or systematic) 
internal correlation between data. Overall, the use of the Kurtosis test only may not lead 
to a conclusive assessment of Gaussianity. 

For this purpose we adopt the Shapiro-Wilk (SW) tests, which is suitable to test the 
departure of a data sample from normality. SW tests whether a data sample {xi, . . . , x n } of 
n observations comes from a normally random distributed population. The so-called Wsw- 
statistic is the ratio of two estimates of the variance of a normal distribution based on the 
data sample. In formula: 

Wsw ~ ELiK-^) 2 ' (20) 

where X(h) are the ordered sample values {xq^ is the h-th. smallest value) and are weights 
given by [46]: 

(";■ ■■■Q>n) _ — •• 



(m T \/- 1 \/- 1 m)2 
12 



with m T the expected values of the order statistics of random variables sampled from the 
standard normal distribution, and V is the covariance matrix of the order statistics. From 
a mere statistical point of view, Wsw is an approximation of the straightness of the normal 
quantile-quantile probability plot, that is a graphical technique for determining if two data 
sets come from populations with a common distribution. Notice that W$w e [0, 1]. One 
rejects the null hypothesis of normality within a significance interval of 0.05, if p-Wsw < 0.05, 
where the p-Wsw is the p- value of Wsw i-e. the probability of obtaining a result at least as 
extreme as the one that was actually observed, given that the Gaussian hypothesis is true. 

The two tests verify two complementary aspects. Even if the SW one is considered 
a faithful Gaussianity test it can fail either for a non-Gaussian or for non truly random 
distributions. Once the Gaussianity of the data is proved, by means of the Kurtosis excess 
7, the SW test is used as a test for the randomness of the data ensemble. 

We have applied the above statistical analysis to our homodyne data distribution divided 
into 104 discrete phase bins (each bin correspond to a ^variation of « 60 mrad). As an 
example of Gaussianity test, in Fig. 2 we show two typical experimental homodyne traces 
for modes b and d (plots on the left) as well as the corresponding p- value of the Shapiro- Wilk 
test (plots on the right). As it is apparent from the plots, the mode b is excited in a thermal 
state, while the mode c is squeezed with quadratures noise reduction, corrected for non-unit 
efficiency, of about 2.5 dB. An analogue behavior has been observed for the other modes. 
Since we have p-value p > 0.05 (the dashed line in the plots) for all the data set we can 
conclude that our data are normally distributed and that the signals arriving at the detector 
are Gaussian states. 

6. Tomographic reconstruction 

As already mentioned in Section 4 our setup is suitable to measure all the field quadrature 
x e = x cos 9 + y sin 9 of any input mode by scanning over the phase of the homodyne local 
oscillator. We exploit this feature twice. On the one hand we use the full homodyne set of 
data to assess Gaussianity of the state and, on the other hand, we may perform full quantum 
homodyne tomography to validate results and increase precision for some specific quantities 
[47]. 

The acquisition of every mode is triggered by the PZT linear ramp: for each value 9, the 
quadrature Xq = x cos 9 + y sin 9 of the homodyne input mode is measured, where x and y are 
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respectively the amplitude and phase field quadratures. Calibration with respect to the shot 
noise is obtained by obscuring the OPO output and acquiring the vacuum quadratures. It is 
worth stressing that experimentally, the acquisition over 2n intervals presents the advantage 
that it does not require sophisticate phase locking set-up to keep 9 constant during the 
acquisition. 

The collection of homodyne data points, normalised to the shot-noise, is then used to 
evaluate the bipartite state properties, included the quadratures, xq, for every 9. The deter- 
mination of the quadratures mean value, as well as of any other relevant quantity, has been 
performed thank to the so called pattern function tomography. This allows reconstructing 
the mean value (O) of an observable O as the statistical average of a suitable kernel func- 
tion R[0] over the ensemble of homodyne data (xi,0i) [48]. By taking into account the 
non-unitary detection efficiency 77, (O) is indeed retrieved as: 

1 N 

(0) = R[d] = -J2 R vP] {Xi-A) 

1=1 

where N is the total number of samples. Every datum (xj, Ot) individually contributes to the 
average, so that the operator mean value is gradually built up, till statistical confidence in 
the sampled quantity is sufficient. Although the method is very general, and can be applied 
to any operator, in the following we will only report the kernels for the quantities we are 
interested in in this paper. For rj > 0.5, the following kernels can be calculated (we omit 
the dependence of R v [0](x; 9) on x and 9) 



I r 

R v [era] = 2x 2 , R v (a' a)' 

2n L 



-x A - 2x z 



R v [x^\ = 2x cos ((f) — 9) 

Rv K] = \{ 1+ ( Ax2 - ^) t 4cos2 x ] 

In principle, a precise knowledge of the (O) would require an infinite number of measurements 
on equally prepared states. However, in real experiments the number of data iV is of course 
finite, so requiring an errors estimation. Under the hypotheses of the central limit theorem 
the confidence interval on the tomographic reconstruction is given as: 



5d = -LJI^d] 
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where AR^[0] is the kernel variance, say the average over the tomographic data of the 
quantity R 2 [0] (x,<p) — (O) 2 . For the particular case of a field quadrature, the confidence 
interval is: 



AR V [x e \ (x, 4>f = (Ax 2 ) + X - (n) + - 



where (n) is the mean photon number of the field under scrutiny. 
7. Experimental results 

The first step is the Gaussianity test for the each data set which consists, for each acquisition, 
of a collection of eight homodyne traces: one for the shot-noise (vacuum), one for the 
electronic noise and six corresponding to the six homodyne modes {a, b, c, d, e, /} Then 



the experimental errors. After this the thermal character of a and b, as expected for a below 
threshold OPO, is verified and their mean photon number as well as A and B CM blocks 
are retrieved. Then, modes c, d and e, /, are analysed in view of their squeezed thermal 
nature, with squeezing appearing on the x, yand t, z quadratures respectively. The variance 
of x c (squeezed), y c (anti-squeezed), Xd (anti-squeezed), ya (squeezed), x e , y e , Xf, and y/, 
are finally used to retrieve the CM C block. 

Since modes a and b are both, phase independent, thermal states, the determination of 
their quadrature variances are highly robust against homodyne phase fluctuations. Accord- 
ingly, the error on blocks A and B elements is obtained by propagating the relative tomo- 
graphic error. On the other hand, when dealing with, phase dependent, squeezed states, a 
small uncertainty in setting the LO phase 9 can result in a non negligible indeterminacy on 
the quadrature variance used to reconstruct the relative cr element. As a consequence, when 
evaluating the errors on the elements of the block C, one must take into account the noise 
properties of the involved modes and critically compare the tomographic error with the error 
due to the finite accuracy on 9. a 13 and <t 2 4 are obtained as combinations of squeezed/anti- 
squeezed variances, which, are stationary points of the variance as function of 9, thus they 
are quite insensitive to 9 fluctuations; accordingly the overall tomographic error can be re- 
liably used in this case. On the contrary au and 023 depend on the determination of x e j 
and y e j. These quadrature variances are extremely sensible to phase fluctuations being the 
variance derivative, in 9, maximum for this values. In this case the error correspond to the 
deviation between the variances at xz. and at xn±se (or x_zl and x-K±se) with 59 ~ 20 



we check the consistency of the vacuum 
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(21) 



mrad, corresponding to the experimental phase stability of the homodyne detection. 

Once the full CM is reconstructed the analysis of the bipartite state can start. First, if 
a failure of the uncertainty condition for the minimum symplectic eigenvalue (see Eq. (6)) 
occurs it means that the measurement must be discharged. In this case the reconstructed 
cr does not correspond to a physical state. 

A typical matrix is given by 

/ 1.694 0.000 1.204 -0.02 \ 
0.000 1.694 0.02 -1.232 
1.204 0.02 1.671 0.000 
\-0.02 -1.232 0.000 1.671 / 

It corresponds to an entangled state that satisfies the Duan criterion (/3 D = 0.93) and 
the Simon criterion (gL = 0.46, E(cr) = 0.12) while it it does not show EPR correlations 
(Pepr = 0.65). This fact is not surprising: the state is rather robust (the mean photon 
number of the system is ps2.4) but the correlation, while ensuring the non-separability of 
the state, does not provide EPR-like effect. Such a state would not be useful in quantum 
protocols where EPR is exploited whereas it is sufficiently correlated for enhancing the 
security of CV QKD. For such a state the entropies are given by S(g) = 2.23, 5*(1|2) = 0.734 
and 5(2 1 1) = 0.720 and the quantum mutual information by /(cr) = 0.779. 

A strongly correlated system is the one corresponding to a diffenre data set and described 
by the matrix 



(22) 



/ 2.107 0.000 1.830 -0.1 \ 
0.000 2.107 0.08 -1.573 
1.830 0.08 1.867 0.000 
\ -0.1 -1.573 0.000 1.867 / 
In this case the corresponding state, whose total energy is n tot 2.9, is both entangled and 
EPR correlated (p D = 0.64, <L = 0.23, E(a) = 1.12, and (3 EPR = 0.22). This reflects in a 
higher value for the mutual information /(cr) = 1.633 carried by t he state. Notice that this 
state suffer from non-zero entries on the anti-diagonal elements of the CM. This is due to a 
non-perfect alignment of the non-linear crystal that give raise to a projection of a residual 
component of the field polarized along a onto the orthogonal polarization (say along b), 
thus leading to a mixing among the modes [45]. This effect is the well known polarization 
cross-talk. 



16 



Indeed, in the ideal case the OPO output is in a twin-beam state S(C)|0), S(C) = 
exp-j^a^ — (ab} being the entangling two-mode squeezing operator: the corresponding 
CM has diagonal blocks A, B, C with the two diagonal elements of each block equal in 
absolute value. In realistic OPOs, cavity and crystal losses lead to a mixed state, i.e. to 
an effective thermal contribution. In addition, spurious nonlinear processes, not perfectly 
suppressed by the phase matching, may combine to the down conversion, contributing with 
local squeezings. Finally, due to small misalignments of the nonlinear crystal, a residual 
component of the field polarized along a may project onto the orthogonal polarization (say 
along 6), thus leading to a mixing among the modes [45]. Overall, the state at the out- 
put is expected to be a zero amplitude Gaussian entangled state, whose general form may 
be written as g = U(/3)S(C) LS(fi, 6) T LS^fi, 6) S t (C)Ut( / 5), where T = n <g> r 2 , with 
Tk = (1 + nfc)~ 1 [^fc/(l + ^fc)]" 1 " denotes a two-mode thermal state with average photons 
per mode, LS(£i,£ 2 ) = £(£1) ® ^(£2), S(£ k ) = exp{|(£ fc a t2 — ^a 2 )} denotes local squeezing 
and U(/3) = exp{(3a^b — fiatf} a mixing operator, (, ^ and (3 being complex numbers. 
For our configuration, besides a thermal contribution due to internal and coupling losses, 
we expect a relevant entangling contribution with a small residual local squeezing and, as 
mentioned above, a possible mixing among the modes. 

Given the CM it is also possible to retrieve the corresponding joint photon number dis- 
tribution p(n,m) by using the relation [7]: 

p(n,m)= f d Al d A2 x(Ai, A 2 ) Xn(-Ai) Xm(~A 2 ), (23) 

where x{^ii^2) is the characteristic of the reconstructed two-mode state, that actually 
depends only on cr, and Xn(Afc) denotes the characteristic function of the projector \h)(h\, 
Xn(A) = (n\D(X)\n) = exp{ — ||A| 2 }L n (|A| 2 ), where L n (x) is the n-th Laguerre polynomials. 
In Fig. 3 we report the joint photon number distribution p(n, m) derived from the CM (22) 
and the single-mode photon distributions (either from data or from the single-mode CM) 
for modes b and d (the same modes considered in Fig. 2): as one may expect, the photon 
number distribution b is thermal, whereas the statistics of mode d correctly reproduces the 
even-odd oscillations expected for squeezed thermal states. 

The mutual information /(cr) of Eq. (9) measures the amount of information one can get 
on one of the two subsystems by measuring the other one. In turn, it is a measure of the 
degree of correlation between the two modes. On the other hand, equally entangled states 
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may show different and the difference appears to be dependent on the total number 
of photons. In Tab. 1 we report the mutual information I(cr), the total number of photons 
n tot , the Duan and EPR factors (3d and (3e, and the symplectic eigenvalue eL for different 
acquisitions. All the states are non-separable and not EPR-correlated; they have different 
number of photons and, correspondingly, different quantum mutual information. 

8. Conclusions 

Gaussian states of bipartite continuous variable optical systems are basic tools to imple- 
ment quantum information protocols and their complete characterization, obtained by re- 
constructing the corresponding covariance matrix, is a pillar for the development of quantum 
technology. As a matter of fact, much theoretical attention have been devoted to continuous 
variable systems and to the characterization of Gaussian states via the CM. On the other 
hand, only a few experimental reconstructions of CM have been so far reported due to the 
difficulties connected to this measurement. 

We have developed and demonstrated a reliable and robust approach, based on the use of 
a single homodyne detector, which have been tested on the bipartite states at the output of 
a sub-threshold type-II OPO producing thermal cross-polarized entangled CW frequency 
degenerate beams. The method provides a reliable reconstruction of the covariance matrix 
and allows one to retrieve all the physical information about the state under investigation. 
These include observable quantities, as energy and squeezing, as well as non observable ones 
as purity, entropy and entanglement. Our procedure also includes advanced tests for the 
Gaussianity of the state and, overall, represents a powerful tool to study bipartite Gaussian 
states from the generation stage to the detection one. 
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List of captions 

• Fig. 1. Experimental setup: A type-II 0P0 containing a periodically poled crystal 
(PPKTP) is pumped by the second harmonic of a Nd:YAG laser. At the OPO output, 
a half- wave plate (A/2 out ), a quarter- wave plate (A/4 out ) and a PBS out select the mode 
for homodyning. The resulting electronic signal is acquired via a PC module. 

• Fig. 2. (Left): from top to bottom, two typical experimental homodyne traces of 
modes b and d (similar results are obtained for the other modes). (Right): p- value 
of the Shapiro- Wilk normality test as a function of the bin number (see the text for 
details). Since we have p- value > 0.05 (the dashed line in the plots), we can conclude 
that our data are normally distributed. 9 is the relative phase between the signal and 
the local oscillator. Kurtosis excess 7 for these data is within experimental error. 

• Fig. 3. (Left): Joint photon number distribution p(n,m) for the entangled state of 
modes a and b at the output of the OPO. (Right): single-mode photon distributions 
p(n) for modes b (top right) and d (bottom right). The single- mode distributions of 
mode b is thermal and corresponds to the marginals of p(n,m). The distributions for 
modes d is that of squeezed thermal state. 
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Fig. 1. Experimental setup: A type-II OPO containing a periodically poled crystal (PP- 
KTP) is pumped by the second harmonic of a Nd:YAG laser. At the OPO output, a 
half-wave plate (A/2 out ), a quarter-wave plate (A/4 out ) and a PBS ou t select the mode for 
homodyning. The resulting electronic signal is acquired via a PC module. 
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Fig. 2. (Left): from top to bottom, two typical experimental homodyne traces of modes b 
and d (similar results are obtained for the other modes). (Right): p- value of the Shapiro- 
Wilk normality test as a function of the bin number (see the text for details). Since we have 
p- value > 0.05 (the dashed line in the plots), we can conclude that our data are normally 
distributed. 9 is the relative phase between the signal and the local oscillator. Kurtosis 
excess 7 for these data is within experimental error. 
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Fig. 3. (Left): Joint photon number distribution p(n, m) for the entangled state of modes a 
and b at the output of the OPO. (Right): single-mode photon distributions p(n) for modes 
b (top right) and d (bottom right). The single-mode distributions of mode b is thermal and 
corresponds to the marginals of p(n,m). The distributions for modes d is that of squeezed 
thermal state. 
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1(a) 


n to t 


Pd 


J_ 


Pe 


0.821 


1.421 


1.54 


0.34 


0.33 


0.480 


1.025 


1.71 


0.42 


0.44 


0.366 


0.879 


1.84 


0.45 


0.46 


0.338 


0.562 


1.67 


0.40 


0.35 


0.220 


0.485 


1.86 


0.44 


0.39 



Table 1. The quantum mutual information for acquisition with different average photon 
numbers together with the Duan and EPR factors (3d and Pe , and the symplectic eigenvalue 
cL. All the states are non-separable and not EPR-correlated. They have different number 
of photons and, correspondingly, different quantum mutual information. I(tr) appears to be 
an increasing function of the total energy of the state. 
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